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We review the status of the understanding of single-electron transport (SET) 
devices with respect to their applicability in metrology. Their envisioned role 
as the basis of a high-precision electrical standard is outlined and is discussed 
in the context of other standards. The operation principles of single electron 
transistors, turnstiles and pumps are explained and the fundamental limits of 
these devices are discussed in detail. We describe the various physical mech- 
anisms that influence the device uncertainty and review the analytical and 
numerical methods needed to calculate the intrinsic uncertainty and to opti- 
mise the fabrication and operation parameters. Recent experimental results 
are evaluated and compared with theoretical predictions. Although there are 
discrepancies between theory and experiments, the intrinsic uncertainty is al- 
ready small enough to start preparing for the first SET-based metrological 
applications. 
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1 Introduction 

1.1 General introduction 

In recent years, the development of sub-micron devices has offered the opportunity to 
transport electric charge by the manipulation of individual electrons. Devices have been 
fabricated that can be controlled by an externally applied alternating signal in such a way, 
that an integer number of electron charges is transferred through the device per cycle. 
Once the frequency of the external signal is accurately known, the current through the 
device can, in principle, be calculated with high precision. Such single-electron transport 
(SET) devices have opened the path towards the realization of high-precision current and 
capacitance standard. A high-precision current standard is needed because a standard 
for the unit of electrical current, the ampere, is presently not available in a direct and 
accurate manner, although it is one of the base units of the International System of Units 
(SI). The construction of a high-precision capacitance standard using SET devices offers 
a realistic alternative for the calculable capacitor. In general, there is also a growing need 
in science and industry to measure very low current levels with high precision, driven by 
the ongoing tendency towards miniaturization of electrical circuits. 

The basis of high-precision current and capacitance standards that are presently be- 
ing developed, is formed by sub-micron devices like turnstiles and pumps. In practice 
the observed intrinsic uncertainty of these devices is of the order of 1 part in 10*, al- 
ready approaching the uncertainty levels of Josephson voltage standards and quantum 
Hall resistance standards. Theoretically, SET devices are capable to individually transfer 
electrons with an error per electron that is much smaller. This discrepancy is caused by 
yet unknown mechanisms. 

Because of the extremely difficult and failure-sensitive fabrication process of SET de- 
vices and because of the ultra-low temperatures needed for their operation, research on 
SET devices is inherently time and money consuming. An efficient feedback between ex- 
perimental results and theoretical understanding is necessary to improve the effectiveness 
of the research activities. For that reason, an overview is desired of today's understanding 
of the properties of SET devices when applied within metrology. 

The focus will be on SET transistors, and turnstile and pump devices. Their working 
principle will be discussed in the first section of this article. A brief introduction to the 
fabrication technique and a short overview of the experimental status are given in Sec. 
|2|. In Sec. ^ and ^ it will be shown that the error mechanisms affecting the accuray of 
SET devices are in general reasonably well understood and that theoretical calculations 
can support the design of new device lay-outs. Nevertheless, in Sec. |5| attention is paid 
to the fact that in many cases the experimentally observed uncertainties are still much 
higher than expected from theory. Finally, in Sec. ^ some alternative devices are briefly 
discussed, such as the recently developed devices that make use of surface acoustic waves. 

Nowadays SET devices are able to accurately generate currents at the level of 1 pA. 
This is already high enough to realize a capacitance standard. Such standard is based on 
charging a capacitor of typically 1 pF with an accurately known number of electrons and 
subsequently measuring the voltage that has developed across the capacitor. However, 
the 1 pA level is still too low to accurately calibrate current sources which operate at 
higher, more convenient levels (microamperes and higher). One of the major tasks within 
SET metrology is to increase the amplitude of the current that can be generated by SET 
devices, while maintaining a very low uncertainty. This article is meant to serve as a useful 
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reference in that process. 



1.2 Historic framework 

The phenomenon of the Coulomb blockade of tunneling provides the possibility to control 
the transport of individual electrons. The Coulomb blockade can be understood from Fig. 
|2| which shows two conducting electrodes coupled to a small island via tunnel junctions. 
Transport of electrons from one electrode to the other can take place by means of tunneling 
through the junctions. However, when the total capacitance Cs of the island is sufficiently 
small, the energy Ec = e^/2CE (with e the electron charge) needed to charge the island 
with an extra electron becomes high enough to prohibit the tunneling of other electrons. 
The Coulomb blockade of tunneling can only be observed at low temperatures k^T <C Eq, 
at which thermally activated transport is suppressed (here is Boltzmann's constant). 

The Coulomb blockade was for the first time experimentally observed in disordered 
granular materials [p. Later, techniques were developed to fabricate — in a controlled 
way — ultra-small tunnel junctions that exhibit the Coulomb blockade effect Nowa- 
days these junctions form the basis of so-called single electron tunneling (SET) devices, 
in which the transport of individual electrons can be controlled by means of externally 
applied signals. Likharev et al. pointed out @ the possibility to use these devices for 
metrological applications such as an electrical current standard, for which the magnitude 
of the generated current is determined by the frequency /set of the externally applied 
signal and the number k of electrons that is transferred per period: 

/set = kefsET- (1) 

The initial suggestion was to pass a current through a single tunnel junction and to 
lock the rate /set/c to an external ac source. The experiments however showed that 
the Coulomb blockade and consequently the sought SET oscillations for a single junction 
are strongly suppressed This is attributed to quantum charge fluctuations in the 
electromagnetic surroundings fs], 0, ^ . 

Special precautions have been taken in several experiments to suppress quantum fluc- 



tuations by attaching highly resistive leads to the junctions |9|, |Tl|]. Effective shielding 
from external quantum fluctuations can also be achieved by placing several tunnel junc- 
tions in series. In this way SET oscillations were experimentally observed in the form of 
current plateaux in the current-voltage characteristics of tunnel junction arrays irradiated 
by a microwave signal However, the flatness of the plateaux was still poor. 

Devices to manipulate electrons in a more controlled way were developed, starting with 
the so-called single electron turnstile |]l3| . A slightly more complicated device is the single- 
electron pump, in which electrons are moved through a series of tunnel junctions without 
applying a bias voltage [^] . The flrst high-accuracy pump was fabricated at the National 
Institute of Standards and Technology (NIST) using a 5-junction design [|^. Recently, an 
improved 7-junction device was fabricated and studied, and an error per pumped electron 



of 1 part in 10^ was achieved [16, 17, 18 1. This accuracy is still much lower than theory 
predicts on the basis of the present understanding of error mechanisms. Further research 
is therefore still necessary. Nevertheless, the achieved accuracy already corresponds to the 
required level for the construction of a new high accuracy capacitance standard proposed 



inRef. p] 



Superconducting devices have also been proposed as a candidate to relate current to 
frequency pO]. Frequency locking at / = 2ef has been observed in small superconducting 
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Figure 1: The metrological triangle connecting frequency, voltage and current through the 
Josephson effect, the quantum-Hall effect and the single-electron tunneling effect. Here 
i,n and k are integers and Rk = h/e^ and Kj = h/2e. 



Josephson junctions (see and references therein), but the effect was strongly sup- 

pressed by thermal smearing and quantum fluctuations 24 1 . Like in the normal metal 
case, it is necessary to consider arrays of junctions in order to eliminate the quantum 
fluctuations p^]. 

Very recently, a device was developed at Cambridge University that utilizes a surface 
acoustic wave (SAW) to pass electrons through a quantum point contact in the pinch-off 
regime ^] . The electric field induced by the SAW leads to the formation of quantum 
dots which move along the constriction transporting an integer number of electrons per 
period of the SAW. The first results are encouraging, but much improvement is needed to 
achieve an uncertainty comparable to that of a SET pump. 



1.3 Metrological motivation for SET research 

The metrological motivation for research on SET devices lies in their expected applicability 
in current standards and capacitance standards. 

A current standard based on the manipulation of single electrons will be the first 
direct and accurate representation of the ampere. Presently the ampere, being the (only) 
electrical base unit within the SI, is represented by the combination of the Josephson-array 
voltage standard and quantized Hall resistance standard through Ohm's law. This is an 
indirect and cumbersome method and in practice only secondary standards for voltage 
and resistance are combined in this way. 

A SET current standard will not be a realisation of the ampere in the sense that 
it is a reconstruction according to the present SI definition, because the basis for this 
definition is Ampere's law, which expresses the force generated by two currents. Ampere's 
law thus connects the electrical and mechanical units within the SI system. Using force 
measurements to realise the ampere is, however, not very practical. A far more convenient 
way appeared to be the application of Ohm's law and to combine voltage and resistance 
standards. In 1990 it was agreed that the quantum-Hall effect and the Josephson effect 
are representations of the volt and ohm, respectively. The assignments of values of the 
new quantum standards were done by experiments where, e.g., the Josephson voltage was 
measured in terms of mechanical units. The resulting relative uncertainties are of the order 
10~^, which is much larger than the reproducibility of the Josephson voltage standards. 
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which is of the order 10^^. The search towards a SET current standard can therefore also 
be viewed in a broader perspective, namely to bring up arguments for a redefinition of 
the SI system, where the electrical units are defined in terms of quantum standards, and 
mechanical units are being derived from those. 

Presently, it is possible to make an accurate SET current standard only for small 
current values, typically a few picoampere. The uncertainty of the current is dominated 
by the uncertainty in the average number k of electrons that is transferred per cycle of 
the drive frequency /set- For a SET standard that can be accurately linked to secondary 
standards, the current level has to be increased by at least several orders of magnitude. 
However, simply increasing the current by means of increasing /set is not a solution since 
this would rapidly increase the relative uncertainty (see Sec. ^). It is an inherent drawback 
of nowadays current standards that are based on SET devices. 

A capacitance standard based on SET devices does not need currents higher than a 
few picoampere, in contrast to a current standard. For an application in capacitance 
standards, the SET device is used to charge a capacitor with an accurately known number 



of electrons [28, 19 1. By measuring the voltage that consequently develops accross the 
capacitor, the zero- frequency capacitance value can be determined. Again, the uncertainty 
of measurement is dominated by the uncertainty in the average number k of electrons 
that is transferred per cycle of the drive frequency /j. It is expected that capacitance 



measurements can be carried out with relative uncertainties of 10""^ at 1 pF |2£]. 

Besides their potential usefulness as measurement standards for calibration of current 
and capacitance, SET devices can also be exploited for fundamental metrology experi- 
ments, such as an alternative determination of the value of the electron charge. One way 
of doing this is to use Ohm's law by passing a SET current through a quantized Hall 
device while comparing the Hall voltage with a Josephson-array voltage standard. This 
closes the so-called quantum metrological triangle. Fig. first suggested by Likharev and 
Zorin |2|. 

The quantized Hall resistance R and the Josephson voltage V are given by 

V = nh/Kj, R = RK/i, (2) 

where n and i are integers, and /j is the drive frequency of the Josephson-array voltage 
standard. The constants of proportionality Kj and Rk are called the Josephson con- 
stant and the Von Klitzing constant, respectively. Practical values for Kj and Rk have 
been recommended to facilitate the international comparison of voltage and resistance 



measurements |30[]. The recommended values are denoted by i^j-go and R^^go, respec- 
tively. By combining Eqs. (||) and (^) and inserting values for Kj and Rk , a value for 
the electron charge can be calculated from 



V/SET/ V i J \RkKj 

Of course, this value has to be consistent with earlier determinations of e. Based on the 
physical understanding of the Josephson effect and the quantum Hall effect, the constants 
Kj and Rk are predicted to be combinations of e and h, 



2e „ h 



= —, Rk = -. (4) 



This means that, in addition, the value for /j//set y^nk/i is expected to be exactly equal 
to 2. It should be stressed that a comparison of the so measured value of e with the 
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already known value of e does not yield any additional information about the values of 
the fundamental constants, or the correctness of the underlying theory. It merely provides 
a usefull consistency check, unfortunately without providing a clue with respect to the 
source of the possible inconsistencies. 

In a more direct way, the electron charge can be measured by charging a capacitor 
Cs with an accurately known number m of electrons, using a SET current source. When 
the value of the capacitor is accurately known, through traceability towards the calculable 
capacitor, and by measuring the voltage across the capacitor using a Josephson standard, 
the value of e can be calculated from 



Again, this value must be consistent with the generally accepted value for the electron 
charge. Alternatively this method of charging a capacitor can be used for a determination 



of the fine structure constant a [1£], which is defined by 

(6) 

where fiQ is the permeability of vacuum and c is the velocity of light. Combining Eqs. 
and the value for a is given by 

"=f^^)xf^)xAC. (7) 



V 2 y \2m 

Such a determination of a would be helpful in the field of fundamental constants. 



1.4 Basic requirement for the Coulomb blockade of tunneling 

The Coulomb blockade of tunneling can only be observed if the number of electrons on 
each island is well defined. First of all, this implies that the energy needed to place an 
additional charge on an island must be large compared to the energy related to thermal 
fluctuations: 

Ec = eV2Cs » kBT, (8) 

where Cs is the total capacitance of the island. Secondly, quantum fluctuations must be 
suppressed. In other words, the wave function of an electron should be localized on a 
single island. Therefore, the energy uncertainty associated with the lifetime RtC due to 
tunneling of electrons across the junction (with capacitance C and tunneling resistance 
Rt) must be much smaller than Eq- From this follows: 

Rt > h/e^ ~ 25.8 kn, (9) 

where h is Planck's constant. 

The latter condition suggests that the larger the tunneling resistance the better. This 
is unfortunately not the case, because the speed of operation is limited by the time RtC, 
which therefore should be as short as possible. One thus has to find an optimum choice 
of Rt which compromises between speed and leakage. 
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Figure 2: SET transistor. The top panel shows a diagram of a SET transistor with a 
mctahc island which is connected to the external electrodes through two tunnel junctions 
(represented by boxes). The capacitances of the tunnel junctions axe Cl and Cr. The 
lower left panel shows the calculated current-voltage characteristic at low temperature 
(T = Ec/lO) and Vg = 0. For a bias voltage above threshold, Vtr = e/Cs, the current 
can flow, whereas for smaller bias voltage, the current is blocked. The number of electrons 
on the island can be tuned by the gate voltage Vg- The properties of the transistor change 
periodically with the gate voltage Vg- The period corresponds to addition of one electron 
to the island. This is shown in the lower right panel, where the conductance is plotted 
versus Vg- The current reaches a maximum each time the charge state with N electrons 
has the same energy as the state with N + 1 electrons. 
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1.5 The SET transistor 

1.5.1 Properties of the SET transistor 

To elucidate the effect of the Coulomb blockade on the electric properties of nanometer 
scale devices we consider the so-called single electron transistor whose diagram is shown 
in the upper part of Fig. |2|. The central island can only be occupied by an integer number 
of electrons, provided that the conditions (^), (|9|) are fulfilled. By applying a bias voltage 
(Vl — Vr) and changing the potential Vg of the gate electrode the number of electrons 
on the island can be changed. If Vg is adjusted such that the charging energy, needed 
to add (remove) an electron to (from) the central island, is larger than the applied bias 
voltage and the thermal energy, current through the structure is blocked. If however the 
potentials are selected so that the charging energy is zero, current can fiow. The device 
thus acts like a transistor where a small change in offset charge on the gate capacitor 
changes the state of the device from closed to opened. For comprehensive reviews of the 
Coulomb blockade see e.g. Refs. @,|3|. 

The electrostatics of the system is determined by the external voltages and the charge, 
Q = Ne, on the central island. The part of the total electrostatic energy which depends 
on the number of electrons is 

Ee^=^-^^^, (10) 

where the total capacitance of the island is Ce = {Cl + Cr + Cg) and the induced charge 
on the island is defined as Qex = {VlCl + VrCr + VqCg) ■ The charging energies for 
adding (A"^) or removing (A~) an electron from the island are given by 

A± = E,,{Q) - E,i{Q ±e) = (q - Qe. ± • (11) 

Thus if Qex = {n + ^)e the two charge states Q = ne and Q = {n + l)e have the same 
energy. As a result, current can flow freely because the system can alternate between the 
two degenerate states. Away from this degeneracy point the current is blocked as long as 
e{VL - Vr) < eV2C72 and A;bT < e^/2Cj:. 

The junctions in a SET transistor normally consist of two layers of metal (usually Al) 
separated by thin insulating barrier ( Al oxide) . The overlap of the initial and final electron 
states \i) and |/) on the two sides of the barrier is described by the tunneling Hamiltonian, 

y = T.Tu\f){i\ : (12) 

where Tj^j is the transition probability from state \i) to state |/). Fermi's golden rule 
allows one to evaluate the rate Tif of tunneling from state \i) to state |/), 



27r 



2 



f V i) 6{Ef-Ei) , (13) 



Using these formulae, we can find the electron tunneling rates through the left and right 
junctions in Fig. |2|, 

1 A± 

Tl.r (g ^ Q ± e) = ^- (14) 

where Rt{l,r) = ^/2'7re^|ri^ijpA^|. are the nominal tunneling resistances of the junctions 



and Np is the density of states at the Fermi level. From Eq. (14) we see that tunneling 
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is exponentially suppressed if the electrostatic energy after the tunnel event is larger than 
before (i.e. < 0). 

It is important to realize that the tunneling rate, given by Eq. (0), relies on the 
assumption that electrostatic equilibrium is restored instantaneously after a tunneling 
event. In terms of time scales this is true if the characteristic time for redistribution of the 
charge, which is the inverse plasma frequency, is much shorter than the "tunneling time" 
h/A. This assumption is sometimes called the capacitance model (which corresponds to 
the "global rule" in Ref. |^T|) and it is certainly a valid approximation for good metal 
conductors. 

The dynamic properties of a SET transistor are determined by the tunneling rates 
dU), and can be described by the master equation 

j^P{Q, t)=J2 [PiQ ± e, t)r(Q ±e^Q)- P{Q, t)r(Q ^ Q ± e)] , (15) 

where F = E^ + F/j is the sum of the tunneling rates through both junctions and P{Q) is 
the probability of the charge state Q. In a steady state {dP{Q, t)/dt = 0) there is a balance 
between the number of electrons entering and leaving the island so that the probability of 
a transition from the state Q to the state Q + e is equal to the probability of a transition 
in the opposite direction. Thus it follows that 

PiQ)T{Q -^Q + e) = P{Q + e)F(Q + e ^ Q) . (16) 

This equation yields P{Q) and once that is found the current is given by 

/ = eJ^PiQ) [^l{Q -^Q + e)- Tl{Q -^Q-e)]. (17) 
Q 

In the lower part of Fig. |2| we illustrate the behavior of the single electron transistor by 
plotting the characteristics based on the solution of Eq. (|l6|). The source-drain conduc- 
tance depends strongly on the additional charge on the island induced by the gate voltage. 
The distance between two conductance minima corresponds to the addition of one elec- 
tron charge to the island. Clearly, a SET transistor is a very sensitive detector of electric 
charge. 



1.5.2 Applications of the SET transistor 

Since the source-drain conductance of a SET transistor is strongly dependent on the 
induced charge on the island it is possible to use the SET transistor as a very sensitive 



detector of changes in gate voltage |32, This can be done by applying a suitable bias 
current and by measuring the dependence of the source-drain voltage V^d as a function 
of the gate voltage Vq- After optimisation of the values of the junction capacitances and 
gate capacitance a voltage gain dYso/dVc close to one can be obtained. This electrometer 
operation of the SET transistor however only works at low frequencies due to inevitable 
stray capacitances. The low-frequency properties are strongly influenced by noise sources 
inside and close to the tunnel junctions. These sources give rise to a 1// - like background 
charge noise on the gate electrode with a magnitude of typically 10~^ e/\/Hz at 10 Hz. 



This will be discussed in Sec. 3.5. 



The advantage of a SET electrometer over conventional electrometers is their very 
low leakage. If the role of stray capacitances can be minimized, SET electrometers are 
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especially useful when applied as a null detector in the experiment of charging a cryogenic 
capacitor using a SET electron pump |jl^. This was first shown theoretically in Ref. 
1 34 1 and has recently been confirmed experimentally [^]. Another application of a SET 



transistor is as a non-invasive probe of the local chemical potential in two-dimensional 



electron gas systems [p8l 3S]. 



The electrometer operation of the single-electron transistor however only works at low 
frequencies due to the parasitic lead capacitances which limits speed. For higher operation 
speed an on-chip amplifier is necessary. However recently a new method was developed 
to avoid this drawback. Schoelkopf et al. [^] have demonstrated a combination of a SET 
transistor and a radio frequency resonant circuit, a so-called RF-SET device. The SET 
transistor was used to modify the damping of the resonant circuit (a similar principle is 
used in RF SQUIDs). This RF-SET device has a charge sensitivity of the order 10~^ 
e/\/Hz at 1 MHz, which no doubt will open for new areas of application of single-electron 
electrometers, see also Sec. 



Pekola et al. showed that a SET transistor can be utilized as a thermometer [41, 42]. 
The idea is to use the thermal smearing of the Coulomb gap in the current-voltage char- 
acteristics as a measure of temperature. In contrast to other single-electron devices, such 
a thermometer can operate at temperatures somewhat larger than the temperature corre- 
sponding to the Coulomb charging energy. The advantage of the SET based thermometer 
is its insensitivity magnetic fields, unlike other cryogenic thermometers. In fact this ther- 
mometer measures directly the ratio of temperature to voltage. 

Recent technological advances in fabrication have resulted in a SET transistor po- 



sitioned at the end of a sharp glass tip |35]. Such a single-electron scanning electrome- 



ter(SETSE) is capable of mapping static electric fields and charges with a spatial resolution 
of 100 nm and a charge sensitivity of a small fraction of an electron charge. The SETSE 
has been used to image and measure depleted regions, local capacitance, band bending, 
and contact potentials on the surface of semiconductor samples. Another application of a 
SET transistor is as a non-invasive probe of the local chemical potential in two-dimensional 



electron gas systems [38, 3£] 



1.6 Single-electron pumps and turnstiles 

Single-electron pumps and turnstiles are devices that generate a dc current satisfying the 
fundamental relation (|^) with high accuracy. They consist of a series of tunnel junctions 
in which the electric potential of the islands are manipulated such that only one electron is 
transferred through the device per cycle of an externally applied signal. The general lay- 
out for both devices is shown in Fig. |3|. The islands with charges Qi, . . . , Qn {Qi = Nie) 
are separated by tunnel junctions with capacitances Ci,...,Cn+i. The electric poten- 
tial of each island can be controlled by the gate voltages Vqi , • • • , Vcn through the gate 
capacitances Cci, • • ■ , Ccn- 

The dynamics of such a device can be described by a master equation for the proba- 



bilities -P(Q, t) of the charge states Q = {Qi, ■ ■ ■ , Qn) of the system [|^, 44, 45[, which is 



a generalisation of Eq. (15) for the non-stationary case 



^p(Q, t) = j2 [m', t)m' - Q) - ^'(Q, t)m - qo] ■ (is) 

Q' 

The first term comes from tunneling into the charge state Q with a rate r(Q' — > Q) 
whereas the second term is due to tunneling out of this state with a rate r(Q — > Q'). The 
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Figure 3: Principle of pump. A series of islands i = 1, . . . ,n, with charges Qi, are connected 
by tunnel junctions and coupled through capacitances Cci to gate voltages. These gate 
voltages control the charges on the islands. By suitable periodic variations of the gate 
voltages a well-defined number of electrons can be transfered through the device. The 
turnstile has a finite voltage drop across the array, and typically only one time-dependent 
gate voltage is applied to the central island. The pump usually operates without a dc bias. 



tunneling rates T are time dependent if is time dependent, see Eq. Numerical 
simulations based on this master equation are reviewed in Sec. ^ 

Notice that we described the single-electron tunneling process in terms of the orthodox 
rate equation (^), which assumes thermalization of the system between each tunneling 
event. If the system does not have time to equilibrate between tunneling events the non- 
equilibrium distribution of the electron system has to taken into account. This has been 
considered in Refs. [36, 37]. 



1.6.1 The electron pump 

The simplest electron pump contains two islands (n = 2) and three tunnel junctions 
(see Fig. |3|). By changing the gate potentials Vgi and Vg2 the charge configuration 
Q = {Qi,Q2) can be controlled. Pump operation starts by increasing the potential of 
the island " 1" such that an electron tunnels onto this island from the left electrode. Then 
the potential of the island " 1" is lowered again, while simultaneously increasing the poten- 
tial of the island "2". This results in tunneling of an electron from the island "1" to the 
island "2". Finally, by lowering the potentials of both islands, an electron from the island 
"2" will be forced to tunnel to the right electrode. The device has then returned to its 
initial charge state and one electron has been pumped from left to right. Typical cycling 
rates are in the range of a few megahertz, i.e. RF frequencies, corresponding to currents 
of the order of a picoampere. For optimal performance the bias voltage {Vl — Vr) must 
be near zero. 

One can illustrate the pump operation by means of the diagram in Fig. ^, which 
displays the regions of minimum Coulomb energy in the {Qgii Qg2) plane for the different 
charge states {Qi-,Q2)i where Qci = CgiVgi^ = li2). The Coulomb energy -E-ei is given 
by 



{Qi-QgiY + {Q2-Qg2Y + {Qi-Qgi){Q2-Qg2) ■ (19) 



If two harmonic signals, shifted in phase by 7r/2, are applied to the gates the system 
will follow a circular path in the {Qgij Qg2) plane. The type of pump operation described 
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above corresponds to the situation where the circle encloses one of the triple points {viz. 
the transfer between the states (0,0), (1,0), and (0,1), as shown in Fig. Q). The path 
in the {Qgi-,Qg2) plane can be chosen in many ways and the optimum choice will be 
discussed in section Sec. 4.3. 



1.6.2 The turnstile 

The simplest electron turnstile contains three islands and four tunnel junctions (upper part 
of Fig. |5[). The turnstile differs from the electron pump by the fact that the external drive 
frequency needs to be applied only to one gate electrode coupled to the middle island. On 
the other hand, a non-zero bias voltage is required to define the direction of the current. 
A single drive signal is an advantage compared to the more complex pump operation. 
The finite bias is a disadvantage since it causes energy dissipation in the device. At the 
very low temperatures, required for the proper operation of SET devices, this can be a 
major problem since the transfer of heat from the metal islands to the substrate is poor 
which causes heating of the electron gas. This makes the turnstile less favourable for high 
precision usage. 

The operation of the turnstile is illustrated in the lower part of Fig. |5|, which displays 
regions of minimum Coulomb energy in the iVcV) plane for different numbers N of 
additional electrons on the middle island. For a finite value of the bias voltage V = 
{Vl — and Vq = 0, the turnstile is in the N = region. The turnstile cycle starts 
by increasing the potential of the gate electrode. When the system leaves the = 
region, an electron tunnels from the left electrode, via the two junctions on the left, onto 
the middle island. This brings the system to the A = 1 state. By subsequently lowering 
the gate potential, the system returns to the A = state as an electron tunnels from the 
middle electrode, via the two junctions on the right, to the right electrode. In this way a 
single electron charge is guided through the turnstile for every cycle of the gate voltage. 
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Figure 5: Schematic layout of the single-electron turnstile device. The device has four 
small tunnel junctions and a gate capacitance. In the lower part of the figure, the stability 
diagram in the (Vg, V) plane is shown. Each square represents a stable configuration with 
N electrons on the middle island. The ratio of the tunnel junction capacitance and the 
gate capacitance has the optimum value: C/Cq = 2. One turn around the operation 
trajectory depicted in the lower part, transfers one electron through the device in the 
direction given by the bias voltage. 
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Figure 6: Illustration of the double angle evaporation technique. Starting with a sub- 
strate covered with resist (A), a pattern is written using electron beam lithograpy. After 
development (B), the structure is etched in several steps, such that a free-hanging bridge 
is created (C). A metal layer is then evaporated (D) from angle 1, and after oxidation 
and evaporation of the second metal layer from angle 2, the material used to form the 
mask is removed and the tunnel junction is created (E). Note that due to the double angle 
evaporation two copies of each feature are produced, as can be seen in Fig. 0. 



2 Experimental status of pumps and turnstiles 
2.1 Fabrication technique 



In order to actually realise the devices that were discussed in Sees. |l.5| and |l.6| , the con- 
ditions set in Eqs. and (^) must be met experimentally. The smallest structures, 
with known geometry and well controlled tunneling resistances, that can be fabricated 
using nanoscale technology yield capacitances of the order of parts of a femtofarad. This 
means that the operating temperatures must be significantly lower than 1 K, which can 
be obtained by the use of a dilution refrigerator. 

The commonly used fabrication technique for SET devices is the so-called double angle 
evaporation technique. A resist layer is applied onto some substrate. Using electron beam 
lithography a pattern is written in the resist layer. After development of the resist the 
resulting nanometer pattern is used as a mask for the evaporation of (usually) Al. The 
evaporation is performed under two angles with an intermediate oxidation step. The 
tunnel barriers are formed by the oxidization of the first aluminum layer. The technique 
is clarified in Fig. ^. Nowadays, tunnel junctions with typical dimensions of 50 x 50 
nm^ can be made more or less routinely. Using this technique, different types of metals 
can be applied to make a varyity of junction types: normal-insulating-normal junctions 
(NIN, e.g., Al/Al oxide/Al junctions exposed to a small magnetic field), superconducting- 
insulating-superconducting junctions (SIS) and superconducting-normal-superconducting 
junctions (SNS). Transport through the latter two types of junctions takes place by means 
of Cooper pairs. 
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Figure 7: Atomic force microscop image of a part of a seven junction single-electron pump 
|18|. The shadow evaporation explained in Fig. ^ results in a structure where each feature 
appears twice (e.g. the two parallel leads leading to the capacitances). The size of the 
island is approximately 700 /xm, and the heigth of each layer is approximately 50 nm 
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Figure 8: Current-voltage characteristics for turnstiles by Geerligs et al. taken from 
Ref. [13|. The left panel shows experimental data. The right panel is a comparison 



with numerical solutions of the master-equation. The lower curve is without applied ac 
potential on the gate and the next ones are for increasing ac power. Dotted lines indicate 
the expected plateau value, I = ef for the 5 MHz applied signal. 



2.2 Uncertainty 



In the early 90's the first turnstiles were made at Delft University of Technology |13, 46, 47] 
and the first electron pumps were fabricated at the Centre d'Etudes Nucleaires de Saclay 

Measurements on a four-junction turnstile is shown in the left part of Fig. |8[ In 
the right part of the figure, calculated I — V curves are plotted. The calculations are based 
only on classical tunneling probabilities (see Eq. ([l^)), not taking into account cotunneling 
(see Sec. ^). The value of the generated current when operating at the plateau agreed 
with I = ef within the measurement uncertainty of 0.3 %. Despite the first promising 
results, no determination of the turnstile current with a higher accuracy has been reported 
since 1990, even though the group in Delft University of Technology later continued the 



work on turnstiles |49| using more advanced fabrication techniques. 

In 1994 a 5-junction electron pump was fabricated at NIST with an accuracy, expressed 



in terms of the error per pumped electron, of 0.5 part in 10^ |15]. The equivalent diagram 
of this pump is shown in Fig. |(a). The sequence of triangular-shaped voltage pulses, 
optimized according to Ref. |5C|, was applied to the four gate electrodes, shown in Fig. 
|9|(b). The pulses were carefully designed to compensate for cross-capacitances. The elec- 
tron pump was used to transfer an electron to and from the neighboring island, so that 
the average charge transport in time is zero. The charge on the island was monitored by 
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Figure 9: Electrical diagram and operational principle of a 5-j unction pump, taken from 
Ref. [^]; (a) shows the equivalence diagram, (b) the sequence of charge polarizations 
applied to the islands (the traces are offset for clarity), (c) the Coulomb energy versus 
an electron on island island n at time indicated by dashed line in (b). The solid arrow is the 
desired tunneling, while dashed and dotted arrows show unwanted thermal or cotunneling 
processes, respectively. 



a SET transitor, which was capacitively coupled to the island. In this way the error rate 
of the pump was determined. 

The next generation of pumps made at NIST contained 7 junctions [l^, In 
Fig. 1^ an atomic force microscop picture of such a pump is shown. For this type of pump 
the error per pumped electron was reduced to 1 part in 10®. This intrinsic error rate is 



already low enough to try to realize a capacitance standard as discussed in Sec. 1.3 



3 Error mechanisms 

There are a number of error mechanisms which deteriorate the high-accuracy operation of 
single-electron devices. Some of these errors are due to classical effects (cycle missing and 
thermally activated tunneling) , whereas others are due to the quantum mechanical nature 
of the electron transport (co-tunneling and photon assisted tunneling), which becomes 
increasingly important at low temperatures. In the next subsections we explain the physics 
of these error events. We will not review the calculation of the combined error rates. For 
those results, we refer to the original articles, e.g. the table in Ref. [^] and the numerical 
works reviewed in Sec. |^. 
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3.1 Cycle missing 

When a single electron device is driven by signals of sufficiently low frequency, the dis- 
tribution of electrons on the islands will at any time correspond to the state with the 
lowest energy. This is called the adiabatic regime. At higher frequencies the electrons are 
not given enough time to tunnel to the state with lower energy and the device operates 
non-adiabatically. In other words, in this regime a so-called cycle missing may occur. For 
example, in the electron pump a tunneling event may be missed when we circulate around 
the triple points as shown in Fig. ^ and the boundaries between the different charge states 
are crossed too rapidly. This means that on the average less than one electron is trans- 
ferred per cycle, which leads to deviation from the ideal behaviour as described by Eq. 

The analysis of a 3-junction electron pump with harmonic drive has been carried out 



in Refs. |51, q^]. It was found that reliable operation is only possible for frequencies much 
smaller than 1/{RtC), where C is the junction capacitance. For an A^-junction pump, 
the error can be estimated by calculating the probability that an electron does not tunnel 
during the time where one of the gate charges, say Qi, raises from to e. This results in 
a error rate per cycle of |50| 



where / is the operation frequency. From these expression it is clear that it is the RtC- 
time that sets the limit for the operation speed. A small RtC time is desirable for high 
speed performance, but, on the other hand, a small tunnel resistance Rt will lead to large 
quantum fluctuations and unwanted tunneling events such as cotunneling, described in 



Sec. 3.2 



In order to minimize the cycle missing error it is important to optimize the waveform 
of the applied signals. In this way the time spent in a state where the expected tunnel- 
ing is energetically favorable, can be maximized. For the turnstile a square-type wave 
form rather than a sinusoidal wave form is preferable. For the pump a triangular signal 



with optimized amplitudes is used in metrological experiments [15, 16, 17, |l^]. Further 



optimization of the waveform using numerical algorithms is reviewed in Sec. 4.5. 



3.2 Thermally activated errors 

Thermally activated tunneling is the main source of errors at high temperatures. In Fig. 
|9|, a series of thermally enhanced tunneling events is indicated by the dashed arrow. From 
Eq. (14) it follows that the tunneling rate is suppressed exponentially if the energy of the 



final state is much larger than the energy of the initial state, 

lA^I 

Tthermal ~ ^^5- exp (- |A=^| /fc^T) . (21) 

One can thus calculate the probability of a sequence of tunneling the events. By taking 
into account the most important sequences, Jensen and Martinis |5C] found approximate 



formulas for the error rate per cycle due to thermally activated processes 



ethermai ~ cexp ( -d-^ ) , (22) 
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Figure 10: Cotunneling event in a SET transistor. Two electrons tunnel in a coherent 
way such that the charge on the island is the same in the initial and the final states. 
Therefore the process overcomes the Coulomb energy barrier, but is suppressed due to 
limited number of available final states in phase space. Furthermore, the probability of 
the two electron process is proportional to the product of the two tunneling conductances. 



where the parameters c and d depend on the number of junctions and the bias voltage. This 
expression was confirmed by numerical simulations. For the parameters of the accurate 
pump described in Sec. |2.2| , the thermal errors dominate for temperatures larger than 100 
mK, but are strongly suppressed for lower temperatures. 

If the system does not have to relax between tunnel events the electron gas heats up 
and consequently the electron temperature is higher than the lattice temperature. The 
elevated temperature is a result of a balance between input power and the relaxation due 



to electron-phonon coupling [|37|, 36 1. The cooling power is decreasing at low temperares 
as a power law approximately as T,VT^, where V is the volume and S is a dissipation 
constant. On the other hand, the input power is proportional to the driving freqeuncy and 
the charging energy, fEc- Using this line of reasoning an estimated electron temperature 



of 80 mK was found in Ref. |37| using typical parameter for metalic single electron pump 
devices, / = 10 MHz and a base temperature of 10 mK. From this consideration it is 
evident that there is a limit to the lowest operation temperature for single electron devices. 

The understanding of the thermal errors due to tunneling is on firm ground. However, 
thermal effects may also be important for processes which involve higher order tunneling 
processes such as thermally activated cotunneling. This is explained in the next section. 

3.3 Cotunneling 



In the master-equation described in Sec. 1.6, we have sofar included only the processes 



where the electron tunneling through a single junction is involved. However, higher order 
processes are also possible One such process is a cotunneling event where two or more 
electrons tunnel simultaneously through several junctions in a quantum-coherent manner. 



In the case of the SET transistor, explained in Sec. 1.5 a two electron cotunneling 
process can transfer an electronic charge from left lead to the right lead such that the 
number of electrons on the island is the same in the initial and the final state. This 



process is illustrated in Fig. 10. There are two possible intermediate states, one where an 
electron jumps through the right barrier first leaving an excess hole on the island and one 
where an electron jumps through the left barrier first. 

There are two fundamentally different types of cotunneling processes: if the state of the 
electron that enters the island is identical to the state of the electron that leaves the island 
the process is called an elastic cotunneling process whereas if the two states are different 
it is an inelastic process. The latter is generally dominant for the metallic systems where 
the island has a large number of electronic states 
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In Section 1.5 we derived the tunneling rate using Fermi's golden rule. The cotunnel- 
ing rates can be derived in a similar fashion except that now we must include processes 
of higher order in the tunneling matrix elements which connect electronic states across 
the tunnel junctions. The effective matrix elements for the inelastic cotunneling process 
depicted in Fig. |l^ can be found using quantum mechanical perturbation theory. Since it 
is a quantum mechanical process energy need not be conserved in the intermediate state, 
i.e., the state immediately after the first electron has tunneled. Furthermore, the differ- 
ent possible intermediate states add coherently and the effective matrix element T (which 
should be entered in golden rule formula (^)) between initial and final state becomes 

!/m^) - 7— ^-gT + — (23) 



efc/ - E 



The electrostatic energy differences between the initial and the two possible intermediate 
states denoted by and E^ correspond to adding or removing an electron from the 
island, respectively. In order to find the rate for transfer of an electronic charge from left 
to right one must sum over all k,k' ,p,p' , and furthermore include appropriate occupation 



probabilities. At zero temperature these integrals can be performed |55|. 

At low temperature and small bias ksT^eV <^ E^ one can approximate the inter- 



mediate state energies in denominators in Eq. ^ by the electrostatic energy differences, 
because energy of the electron-hole pairs is limited by temperature and bias voltage. In 



this case one finds for the cotunneling rate through the device from left to right |56| 



cotun 



GLGRh M l 
127re3 \E+ E' 



{2^kBTf + {eVf 



V 



1 -exp(-ey/A;Br)' 



(24) 



where V is voltage drop across the device. 



The approximate formulae Eq. (24) and the zero temperature result has been tested 
in several experiments [|^, |5^, 59, 60 1 and very good agreement was found. However, 
outside the range of validity for the approximations, i.e. in the cross-over between the 



cotunneling dominated and sequential tunneling dominated regimes the theory [^Sj 56] 
gives unphysical divergences in the cotunneling rate and needs to be modified. 

To repair the possible divergence in ( p3| ) several authors |55, 62, 63, ^] have 
included a life time broadening due to tunneling in and out of the island . The life 
time broadening results in a imaginary term in the denominators of Eq. ( p3| ) and in this 
way a gradual cross-over between the cotunneling and the sequential tunneling limits is 
obtained. Comparison with experimental data has shown that good agreement is obtained 
for not too small tunnel resistances In passing we note that recent work by 

Konig, Schoeller and Schon ||6^ has shown that a thorough second order perturbation 
theory, which does not suffer the divergence problem, gives very good agreement with 



experiment even for rather small tunnel resistances |67]. Unfortunately, these approaches 
are complicated and it is not clear at this point whether generalization to the cotunneling 
through many junctions can be made sufficiently effectively for practical numerical analysis 
of multi-junction devices. In order to have a simple way of evaluating the cotunneling 
contributions, Jensen and Martinis |^0| have proposed the following approximation scheme. 
They evaluated the integration over intermediate states by distributing energy difference 
between initial and final states equally among the electron- hole pairs. This means that 
for the matrix element in Eq. (^), Ck — Cp = Cpi — e^i = {Ei — Ef)/2. Fonseca et al. 



improved this approximation by shifting the lowest order tunneling rate in Eq. (14), to 
lower energies. By doing this these authors obtained a smaller difference between the 
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approximate tunneling rate and the one obtained from the more correct higher order 
approximations citepasq93,koro92,naza93b,aver94,flen97,lafa93. 

The expression for several consecutive cotunneling events was derived by Averin and 



Odintsov |55]. The Jensen and Martinis and the Fonseca et al. approximations are easily 
generalized for the multiple cotunneling case and they are used in the numerical simulations 
of pumps and turnstiles, reviewed in Sec. 13. We give here the final expression for the a 




multiple cotunneling events of order n |5C], because it is the dominant source of error at 
the lowest temperatures: 

r£L = ^ ( TT S ) S'F^iAE^), (25) 

where Gi is the conductance of tunnel junction number i, AE is the energy difference 
between final and initial states, and 

S= E CliwW^XWj?! (26) 



permutations 



AEk - AE/n 



where the sum is over all permutations of cotunneling sequences going from initial to final 
state, and finally 

^"^^^^ W^^- exp(AABT)-l- 

From this expression we see that the higher order cotunneling events are suppressed by 
the ratio of the junction conductances to the quantum conductance e'^/h to the power n. 
For typical junctions with resistances in the 100 kfi range the n order cotunneling process 
contains a factor (0.04)" due to the first product of Eq. (|25|). However, in contrast to 
the thermally activated tunneling it is not exponentially suppressed and the cotunneling 
leakage events will thus dominate at low temperatures. 

One way of suppressing cotunneling is to increase the number of Coulomb islands. 



Another method which was suggested in Ref. [39| is to increase the resistance of the leads 
that couples to the device. The authors of ref. [ 39 1 showed that for a two junction circuit the 
cotunneling rate the cubic voltage dependence in Eq. (|2^) is replaced by Fcotun ~ y3+22;^ 
where z = Re^ /h is the normalized resistance of the leads. This idea is now being pursued 
experimentally |70|] 

Finally, we mention the so-called elastic cotunneling, which means that energy is not 
dissipated to electron hole pairs in the tunneling process. Consequently, the electron states 



p and p' in Fig. IC are the same. In other words, it is "the same electron" that enters 
and leaves the island. The probability of elastic cotunneling is proportional to the single 
particle level spacing [^]. Thus the process is generally very weak for metal islands. 
However for semiconductor devices this effect indeed gives significant contributions to the 
electron transport. 



3.4 Photon assisted tunneling 

When a SET device is irradiated by an external electromagnetic field, the potentials on 
the electrodes will oscillate. Tunneling events which are otherwise not possible may now 
occur through absorption of photons from the external source. Only small amounts of 
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electromagnetic radiation are needed to introduce considerable errors due to photon as- 



sisted tunneling |71]. Therefore careful attenuation of the high frequency noise is necessary 
by means of shielding, thermally anchoring all leads in the cryostat, and applying copper 
powder filters [72|, lithographically made miniature filter ||7^ or lossy coaxial cables [^]. 



When the device is driven by a periodic signal as in the pump or turnstile, it is possible 
that the drive itself induces photon assisted cotunneling events. This effect was considered 



in Ref. |64], but it is was found that only at a frequency higher than several gigahertz this 
error contribution becomes significant. It can thus be neglected for present days pumps 
and turnstiles operating at typically 5-10 MHz. 

3.5 Background charges 

Noise introduced by a random motion of charges in the dielectric surrounding (the sub- 
strate, the barriers of the tunnel junctions, etc.) of the conducting islands is a serious 
problem in SET devices. This noise leads to fluctuations of the island's polarization. 
Usually, the frequency dependence of the noise spectral density is close to 1//, although 
sometimes it exhibits an almost 1//^ dependence in combination with a telegraph-like 
time dependence E^, 75]. This switching noise results from one or several two- level fluctu- 



ators with typical switching times ranging a fraction of a second up to minutes. For most 
of the aluminium-based SET devices, fabricated on various substrates (Si02, Si, AI2O3) 
the intensity of the background charge noise at 10 Hz is between 10"^ and IQ-^e/V^ 



i|,^|7|,0. 



However, in the devices fabricated with a stacked design the islands are partly |78 



79, ^] or entirely |79, 8C] screened from the substrate by an underneath gate or by an 
underneath counter electrode. Such devices exhibit considerably less noise, down to the 
lowest value measured so far, 2.5 x 10~^ e/\/Hz at 10 Hz |f9| , |80t| . The dominating role of 
substrate in causing the background charge noise is also proven by the mutual correlation 
observed in the cross spectrum of the 1/f noise of two SET electrometers which were 
closely positioned on the same substrate |7^, |8l|. Therefore, a suitably chosen substrate 
and/or an appropriate geometry for SET devices should reduce the background charge 
noise. 

The effect of background charges, as discussed above, is not limited to time dependent 
fluctuations. The dc offset charges |]l5|, |l^ can also arise due to charges trapped inside 
the substrate and the barriers of the tunnel junctions, and due to different work functions 
of dissimilar electrode materials 

4 Numerical analysis of single electron devices 
4.1 Methods and tools for numerical analysis 

Analytical treatment of the single-electron transport is only possible for circuits with few 
tunnel junctions such as the SET transistor, or symmetric devices such as one-dimensional 
tunnel junction arrays. In the last few years, several computer codes have been developed 
to analyze complex single electron circuits numerically |83, 34, |8^, Two main 



approaches have been used. Firstly, Monte Carlo algorithms [|^, gg, g4|, which 
average over various realizations of stochastic SET events. Secondly, algorithms based on 
a numerical solution of the master equation (|l8|), which describes the dynamics of SET 
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circuit in terms of time-dependent probabilities of various charge configurations l5^, 51 



5|, ii,!!, n, g. 



The advantage of Monte Carlo routines is their relative simplicity and robustness. This 
makes them convenient tools for the study of the dynamics of single electron systems, for 
which rare events are not important. Another feature of the Monte Carlo method is the 
trade-off between accuracy and simulation time: one can quickly achieve approximate 
results for very large circuits, which otherwise would not be possible to compute. An 
intrinsic drawback of Monte Carlo algorithms is their inefficiency in the analysis of rare 
events, simply because these rare events might not be encountered during the simulation 
time. Specifically we mention the codes MOSES g and SIMON @ which are able to 
analyze arbitrary circuits consisting of tunnel junctions, capacitances, resistors, and signal 
sources. The latter code is convenient because of its graphical user interface. 

Reliable simulations of SET devices require realistic estimates of the contributions of 
rare events. This can be achieved by starting from the master equation (|T8[), whose solution 
contains complete statistical information about the dynamics of the system including the 
probabilities of the rare events. The problem with the master equation method is the 
necessity to trace the time-evolution of a large number (up to several thousands) of charge 
states. This makes the computation time-consuming. The analysis is alleviated in the 
case of adiabatic devices (like single electron pump) where the electron transfer occurs at 
low energies compared to Ec and the master equation involves a moderate number of low 
energy charge states |5^, |5l|, |9^. In order to use the master equation approach for general, 
non-adiabatic devices (like turnstiles) the development of more sophisticated algorithms 
is required. In the following section, we will briefly describe an effective algorithm of this 
type that is based on the ideas of global analysis and the dynamical choice of the basis of 
states H, 

Finally, we mention the attempt to combine the Monte Carlo approach with solution 
of the master equation |86|. An essential assumption of this method is that the system can 
stay in a rare state for a period much shorter than the characteristic time of the simulation 
(e.g. the period of an external signal in case of pumps). This assumption in fact means 
that the errors due to rare events are analyzed locally in time. In the following section it 
is show that in order to detect errors in the operation of SET devices, the errors should 
be analyzed globally, throughout the whole cycle of the device operation. 



4.2 SENECA: a global approach to rare events 

In this section, we introduce SENECA (Single Electron NanoElectronic Circuit Analyzer), 
a computer algorithm suitable for studies of the dynamics and statistics of single electron 
systems consisting of arbitrary combinations of small tunnel junctions, capacitances, re- 



sistors and voltage sources [35, 53, 54]. The method is based on the numerical solution 
of a linear matrix equation for the vector of probabilities of various electric charge states 
of the system, with iterative refining of the operational set of states. The method is able 
to describe very small deviations from the classical behavior of a system, caused by the 
physical error mechanisms discussed in Sec. ^ . Below, the main underlying ideas and 
characteristics of the algorithm are briefly discussed. 

Global analysis of errors. It must be emphasized that one should distinguish between 
a rare tunneling event and the error in operation of a device (i.e. deviations from Eq. 
(|l|)). Some rare tunneling events bring the system to an unwanted charge state, which 
afterwards evolves classically to a correct operational cycle. These rare events do not 
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affect the proper operation of the device. In order to separate those events from the 
error introducing events, one should analyze the operation globally by solving the master 
equation throughout the whole cycle. 

Clustering of the charge states according to their probability - iterative approach. For 
single electron devices of interest the number M of islands ranges from 5 to 30. If we 
consider the charge configurations of the islands with additional charges and ite, the 
number of states can be roughly estimated as 3*^. Nevertheless, when the temperature is 
low and the tunnel conductances of the junctions are small, most of the possible transitions 
(either those leading to an increase in the electrostatic energy, or those of high cotunneling 
order) have very low rates. As a consequence, a large number of charge states have very 
low probabilities P and may be ignored. Hence, an iterative approach is chosen which 
considers for each iteration n only the limited number of states for which P > P{n). The 
threshold P{n) decreases rapidly with n. 

Dynamical choice of the basis of states. Each iteration corresponds to the solution of 
the master equation at the time interval of the operation cycle. This interval is divided 
into enough time steps, so that the variation of external parameters between steps can be 
neglected. At each time step the algorithm first guesses the set of states whose probabilities 
are expected to be large enough, P > P{n), then solves the master equation for this set, 
and finally filters out the states whose probabilities at the end of the step turned out to 
be low, P < P{n). In this way the minimum possible basis of states is kept throughout 
the computation. 



4.3 Results of simulations 

Pump. First we consider an A'^-junction pump driven by triangular voltage pulses with 



amplitude ae/Cci and offset ue/Cci applied to the gates as shown in Fig. 11(a). The 
gate capacitances are assumed to be small: Cci <C C. If the pump works properly, an 
electron tunnels sequentially through all junctions of the array following the propagation 
of the pulse. The region of correct classical operation (disregarding any rare processes) can 
be evaluated from the condition that tunneling of an electron occurs only in an "active" 
junction biased by a voltage pulse, while the other junctions are closed. By choosing u 
equal to (1 — a)/N one maximizes the region of correct classical operation (Fig. p](b), 

, , (A^-l)(A^-2) r / ^^ 

1^1 < ^ ^ min[a, (2 - a)] , (28) 

with V being the bias voltage in units of e/C 

The closer the parameters a and v are chosen to the center of the region (p8[), the 
smaller the error due to rare events. In particular, the order n of cotunneling processes. 



which, according to Eq. (25) primarily determines the intensity of cotunneling, increases 
from 2 near the inner borders of this region to A" — 1 in its central part (dashed area in 
Fig. 0), 

< ^^^min[a,(2-a)] , (29) 



in which the best performance of the device occurs. In the region (29) the deviation 
AI = I — ef of the current / from its ideal value ef is determined by two competing 
cotunneling processes, shown Fig. |9|. In one process instead of regular tunneling through 
an active junction (corresponding, for example, to the transition (0, 0) (0, 1) in Fig. 
the cotunneling occurs through all the other A^ — 1 junctions in the opposite direction 
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((0,0) —>■ (1,0) —>■ (0,1), (1,0) being the intermediate state, see Sec. ^^ ). In this process 
the electron charge is transferred through the device in the direction opposite to the 
current. The competing cotunneling process starts after the regular tunneling event and 
transfers an electron charge in the direction of the current ((0, 1) — > (1, 0) — > (0, 0) in Fig. 



In the region of optimal operation ( p9D the error AJ in the current can be evaluated 
analytically ||5^, Q . The performance of the device can be characterized by the slope of 
the current platoe, G = d{AI)/dV at 1^ = 0. Being multiplied with the size e/C of the 
optimal region (^) this gives characteristic scale of the current deviation A/* , or accuracy 
of the device, e = Ar/ef = G/fC. 

We consider first the case of a slow drive, ksT > eVo = e^[{N - l)qiRT/NG]^/^ , 
(ji ~ Nf, Rt being the tunnel resistance of the junctions. In this adiabatic regime the 
probabilites to find an electron on the left and right electrodes of an active junction 
correspond to equilibrium conditions at any time. The deviation |A/| of the current due 
to cotunneling increases with temperature and does not depend on driving frequency. The 
accuracy is given by 

Ne^Rr ( irfi 2,, ^>2iv- 



where the factor S (given by Eq. ^ with AE = 0), can be found explicitly in the center 
(a = 1, u = 0) of the optimal region (|29|), 



In the opposite case, fc^T <C eVo the tunneling through an active junction occurs 
non-adiabatically (effectively only a single tunneling event occurs during the gate voltage 
pulse); the deviation |A/| of the current due to cotunneling increases with increasing 



driving frequency and does not depend on temperature |9C, 52|. The accuracy can be 



evaluated as ||9^, 50 1, 



GVq (2A^-3)! \ 7ri?T 



(32) 



Hence, for a given temperature, the optimal driving frequency corresponds to the crossover 
point, where the characteristic energy eVo of electron tunneling is equal to the energy ksT 
of thermal fluctuations. 

So far we have not considered the rare events due to thermally activated classical tun- 
neling whose probability decreases exponentially with temperature, Pth ~ exp(— AE'/Zc^T), 
AE ~ Ec- The analysis [^0[ has shown that thermal activation of classical tunneling can 
be neglected compared to the cotunneling at temperature of the order of or smaller than 
0.025£'c', which is usually the case in experiments. 

The early estimates |5^, |9^ of the abovementioned factors confirmed by more recent 
simulations using SENECA (see Fig. 12) |53, 54] predict that already the 5-j unction pump 
with typical experimental parameters, Rt = 300 kJ7, C = 0.1 fF, / = 10 MHz, T = 
100 mK should have metrological accuracy e ~ 2 • 10^^^. It should be noted that both 
the cotunneling and thermally activated tunneling can be effectively suppressed further 
by increasing the number of the junctions and decreasing the junction size (and/or the 
operation temperature). 
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Figure 11: (a) The shape of the modulation of the gate voltages Va {qi = CGiVa/e) 
in the pump device, (b) The region of correct classical pump operation in (o, v) plane 
(see text). The dashed area is the region where the lowest order of allowed co-tunneling 
transitions is — 1 . Prom Ref . [^] . 

Another way of improvement of the pump operation is the optimization of the driv- 



ing voltage pulse shape [53, 54]. The idea is based on the consideration of the en- 
ergy dependences of the classical tunneling, Eq. (jlj), Td oc AE and cotunneling rates 
Tcot oc AE^^~^, Eq. (p5[). To minimize the ratio Tcot/^d one has to decrease AE during 
the "active" portion of the time interval when the electron transfer between the islands 
occurs. When the probability of the transfer is close to unity, one increases AE to en- 
sure the full transfer of an electron during the "passive" portion of the interval. Being 
optimized by SENECA algorithm, this two-step drive improves the calculated accuracy 
of 5-junction pump (Fig. [T^ ) by a factor of 30 for ksT ~ 100 mK and other parameters 
given above (the improvement enhances at lower temperatures). 

Turnstile. The diagram of correct classical operation of the turnstile with 2A^ junctions 
is shown in Fig. |l3| (see also Fig. ^). Being driven by the gate voltage alternating between 
the two dashed regions in Fig. |l3| the device transfers one electron per cycle of the gate 
voltage. The regions are defined by |9^] 

u + v > {l + c){N -1), u-v <{l + c){N -1) (33) 
u + v < {3 + c)N -1-c, v-u<{l + c){N -1) (34) 

and similar pair of equations with u —u. Here u = {2CgVg/e—l)N, v = {2C+NCg)V/e, 

C = NCg/C. 

Analysis ||9^ shows that the order of cotunneling transitions reaches its maximum 
riQ = 0.8A^ {N ^ 1) for c < 2/3 and the voltages near the threshold of correct classical 



operation {uhigh+v = (l+c)(A^ — 1), ui^w — v = — (l + c)(A^ — 1), see Fig. 13). This happens 
at the optimal voltage Vopt = (1 + c)(no — l)e/[2(2 + c)C]. Careful numerical optimization 
|53, p3] of the accuracy of 6 and 8 -junction turnstiles as function of the gate capacitance 



c and the amplitude Uhigh = —uiow of the gate voltage has shown that the optimal gate 
capacitance corresponds to c ~ 2/3. The calculated uncertainty of 8-junction turnstile 
with optimized parameters was better than 1 part in 10^ for Rt = 300 kO, C = 0.1 fF, 



/ = 10 MHz, and T = 100 mK (see Fig. 0) 
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T(mK) 



Figure 12: Accuracy of the 5-j unction pump driven by optimized step- like waveform (solid 
line), triangular waveform with optimized tunnel resistances (dashed line), and triangular 
waveform with fixed tunnel resistances Rt = 300 kJl (dotted line). The insert shows the 
accuracy as a function of driving frequency at T = 0. From Ref. [^. 
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Figure 13: The diagram of the turnstile operation in the v — u plane, v = (2C + NCG)V/e, 
u = N{^CgVgI e — 1). In the upper dashed region an electron is pulled in the turnstile, in 
the lower dashed region it is pushed out. The central square corresponds to the Coulomb 
blockade state where all classical transitions are suppressed. From Ref. ||90t| . 



Hybrid pump-turnstile. In order to reduce the number of parameters (dc gate voltages, 
amplitudes and phases of RF gate voltages, bias voltage) controlling the operation of a 
pump device, it has been proposed to apply RF signal to every second gate of the pump. 



The calculated uncertainty of a properly biased pump-turnstile with 6 junctions |53, 54] 
was close to the uncertainty of unbiased pump (2 parts in 10^^ for the parameters listed 



above), but deteriorated to 5 parts in 10 for unbiased device, see Fig. 14 



5 Comparison of theory and experiments 

5.1 Discrepancies between measured and predicted errors 

The experiments reviewed in the Sec. ^ have demonstrated an error per pumped electron of 
the order of 1 part in 10^ for pump-like SET devices. Up to very recently, the development 
followed the most obvious path, namely the improvement of the accuracy by increasing 
the number of tunnel junctions (which suppresses cotunneling) and the reduction of their 
size (which enhances the Coulomb blockade). In particular the 7-junction pump has been 
very well characterized |18| with respect to the capacitances of the individual junctions. 



the offset charges at the islands, as well as the electron temperature. 

Substantial efforts have been made ||l5|, |l^ to compare the observed accuracy of the 
pump device with theoretical predictions and to identify unforseen error mechanisms. 
Some details of the experiments are given in Sec. |2.2| . 

For the 5-junction pump |]l^] the error per pumped electron at 40 mK has been ob- 
served to be exponentially dependent on pump frequency for frequencies higher than 5 
MHz. This behaviour is expected when cycle missing is the dominant error mechanism 




Figure 14: Deviation of the current from the quantized value in the unbiased (soHd hne) 
and biased (dotted hne) 5-j unction pump , unbiased (dot-dot-dashed hne) and biased 
(dashed hne) 6-junction hybrid pump (ah driven by a triangular waveform), and 8-junction 
turnstile driven by a step-like waveform (dot-dashed line) at / = 10 MHz as functions of 
temperature. Circles: current deviation at the inflection point Vq, triangles up (down): 
current deviation at Vo(l it 0.1). From Ref. |54|. 



Towards single-electron metrology 



31 



(see 3.1). However, for frequencies smaller than 5 MHz, the observed error per pumped 
electron at 40 mK saturates at a level of 1 part in 10^. This saturation, when ascribed to 
thermally activated processes and cotunneling, occurs at a value that is 10* times larger 
than expected. In order to explain this very large discrepancy Martinis et al. [^] sug- 
gested that the excess error is caused by an additional source of energy, e.g. by leakage of 
thermal noise generated at 4 K or by energy released by charge traps in the substrate. 



For the 7-juction pump |1£], the measured error per pumped electron at 140 mK, when 
the dominant error mechanism is thermally activated tunneling, has been found to be 1 
part in 10^. This is in reasonable agreement with the theoretical estimate of 1 part in 10^. 
However, at 40 mK, when the dominant error mechansm is cotunneling, the measured 
error per electron amounts to 1 part in 10*, which is in sharp contrast to the theoretically 
predicted value of 1 part in lO^''. Since the error per electron is constant and only weakly 
dependent on the wait time between the transfer of individual electrons, the errors are 
concluded to occur predominantly during the active pumping cycle. Therefore, the error 
mechanisms in the pump mode and hold mode are likely to be different. 

Despite the fact that this 7-j unction pump exhibits a much smaller error per electron 
than the 5-junction pump, this improvement is not as drastic as expected from theory, 
especially in the temperature range below 100 mK. 

5.2 Possible sources of additional error mechanisms 

As mentioned in the previous section, additional error sources can exist due to 4 K radia- 
tion or trapped charges in the substrate. Usually, in order to prevent 4 K radiation from 
reaching the device, special microwave filters are placed in series with every electrical lead. 
In addition, the device is mounted in an electromagnetically shielded box. Surprisingly, 
it was shown [jl^ that removal of the filters and the shield does not necessarily affect the 
error rate of the 7-junction pump. This hints at a strong internal low-temperature source 
of spurious tunneling events, which is much stronger than the external noise. It is also 
possible that photons are generated by non-equilibrium charge fiuctuations in the sub- 



strate and/or the tunnel barriers (see Sec. 3A and Ref. |91|). Such an intrinsic photonic 



source cannot be filtered, and may present a fundamental problem for metrological and 
digital application of SET devices. 

Note that it can be difficult to separate the effect of high frequency photons due to 
the charge fluctuations from the effect of low frequency noise due to the drift of the 
offset charges. However, in the experiments reported by Ref. |]l8[ the comparison with 
theory was made when the optimal background charge compensation was stable. Often 
the devices were kept at millikelvin temperatures for several weeks, after which movement 
of background charges was very small and adjustments of the dc voltages on the gate lines 
was only necessary with time intervals of tens of hours. 



6 Alternative designs for single-electron current standards 

The main disadvantage of the discussed current standards based on normal SET junctions 
is their relatively low operation frequencies of a few megahertz, which correspond to a 
current level of only a few picoampere. This limitation is caused by the necessity to have 
large tunnel resistances for the normal junctions Rt ^ Rq = h/e^ in order to fulfill 
the condition of the charge quantization. The limitation in frequency may however be 
circumvented in several ways, as will be discussed in this section. 
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Figure 15: Schematic outline of a SETSAW device, (a) shows the two 2DEGs separated 
by a narrow constriction. Propagating SAWs are indicated by dashed hnes. (b) illustrates 
the potential landscape along the center of the device. In the constriction region a po- 
tential barrier is present. The SAW induces a propagating modulation of the potential. 
In the 2DEG region the modulated potential induced by the SAW is screened and the 
resulting potential is shown in (c) for two different times. In the left part of the figure an 
electron scooped up by a minimum in the potential after which it is moved upward into 
the constriction region (right part of figure). 



6.1 Surface acoustic waves as carriers of single electrons 

The interaction of surface acoustic waves (SAW's) with a two-dimensional electron gas 
(2DEG), present at the interface in a GaAs-AlGaAs heterostructure, can be applied to 
generate an electrical current that is determined by the frequency of the SAW. In Fig. ^ 
the working principle of the so-called SETSAW device is outlined. Due to the piezoelectric 
effect the SAW creates a propagating modulation of the electrostatic potential through 
the 2DEG. Normally this modulation is screened by the electrons in the 2DEG. By using a 
split-gate technique a quasi-one-dimensional channel is formed in the 2DEG. By applying a 
sufficiently negative voltage to the split gate the electrostatic potential in the constriction 
is increased and the electron density can be reduced to zero (pinch-off regime). In this case 
a propagating modulation in the channel is no longer screened. For a sufficiently strong 
amplitude of the SAW the quasi-one-dimensional electron gas in the constriction is split 
into a sequence of small islands, or quantum dots, moving with the velocity of the SAW. 
Recently it has been shown that such a technique can be used to transport an integer 



number of electrons per cycle of the SAW frequency l2q, 27, 93]. The number N of electrons 



per dot determines the dc current / = Nef through the constriction. The experiments 
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have demonstrated the quantization of the dc current at the values corresponding to 
integer N = 1...4. The robustness of the current quantization with respect to a range 
of parameters (gate voltage, bias voltage, amplitude of SAW) implies that the number of 
electrons in each quantum dot is a deterministic rather than stochastic quantity. 

The quantized current can only be observed deep inside the pinch-off regime. The 
concentration n = N/Xsaw of electrons in the pinch-off region, where Xsaw ~ l^ni is the 
wavelength of the SAW, is much smaller than the concentration hq = 4/Xf of electrons in 
the 2DEG, where Ai? ~ 50 nm is the Fermi wavelength. For this reason, the electrons in 
the constriction do not screen the electrostatic potential induced by the SAW. Moreover, 
the kinetic energy of electrons in moving quantum dots is very small, so that the Coulomb 
interaction is the dominant energy scale. The quantum dot that has just been formed 
near the entrance of the constriction can be considered as an electron box connected to 
the 2DEG reservoir. The Coulomb repulsion determines the number of electrons in the 
dot so that the free energy of the system is minimized. The number of electrons in the 
dot can be controlled by the gate voltage and by the amplitude of SAW (the increase of 
the SAW amplitude makes the potential minima deeper which has a similar effect as an 
increase of the gate voltage). 

Due to the substantial difference between the velocity of the SAW, vsaw — 3 x 10^ 
m/s, and the Fermi velocity, vp — 10^ m/s, the formation of new quantum dots occurs 
adiabatically. In contrast to a conventional electron box, the barrier between the quantum 
dot and the reservoir varies in time so that its transparency decreases by many orders of 
magnitude while electrons fill the quantum dot. The dynamics of this process is yet to be 
understood. An important parameter that determines the maximum operation frequency 
seems to be the time dependence of the transparency (or conductance) of the barrier that 
separates the most recently filled dot from the 2DEG. 

The experimental observation of quantized currents at frequencies around 3 GHz sug- 
gests that electron transfer occurs through a low energy barrier, although the accuracy of 
the device may be determined by residual thermal or quantum transport through a high 
energy barrier. Finally we note that another physical mechanism which might negatively 
affect the accuracy - impurity scattering of the electrons in the moving quantum dots - 
becomes less effective with increasing driving frequency ||27[] . 



6.2 Phase locking of coherent Cooper pair tunnehng 

A current standard based on the tunneling of single Cooper pairs may be an interesting 
alternative for the single-electron pump. By using the Coulomb blockade state for Cooper 
pairs, one can take advantage of the coherence which is intrinsic to the superconducting 
state. However, achieving the Coulomb blockade for Cooper pairs is very difficult since the 
effect can be masked by either the onset of a supercurrent or by unwanted quasiparticle 
tunneling. Quasiparticle tunneling is presumably the main cause of the poor performance 
of Cooper pair pumps [|^, |93|| . 

Recently the blockade for Cooper pair tunneling was demonstrated in one-dimensional 
series arrays of small-capacitance tunnel junctions with tunable Josephson coupling |2^. 
In this system the charge is transferred by Cooper pairs solitons which can be considered as 
the dual analog of the vortices in a one-dimensional parallel Josephson array. In particular, 
the Coulomb blockade voltage corresponds to the threshold voltage at which Cooper pair 
solitons enter the array. In the absence of dissipation, the Cooper pair solitons are expected 



to propagate coherently through uniform arrays |94]. A word of caution is needed here. 



Towards single-electron metrology 



34 



In realistic arrays the propagation of Cooper pair solitons might be hindered due to the 
pinning by offset charges. The increase of the Josephson energy effectively reduces this 
pinning. On the other hand, at Ej ~ Ec the onset of superconductivity occurs in the 
array, which destroys the Cooper pair solitons. 

6.3 Electron counting by a high-speed electrometer 

The single-electron current standards discussed so far are based on the active manipulation 
of electrons or Cooper pairs, in other words current generation techniques. A tempting 
alternative is to try to passively count individual charge quanta passing through a voltage 
biased one-dimensional array of tunnel junctions using a high-speed electrometer coupled 
to one of the islands of the array. 

Due to correlations in the single-electron tunneling events, solitons of charge e (or 2e) 
will be running through the array. Mutual repulsion keeps the moving solitons at a certain 
average distance, so that the system of solitons is robust with respect to the random offset 
charges of the islands. Likharev et al. |87| have calculated the time evolution of the 



electric charge at an island of the array. In the frequency domain narrow peaks show up, 
positioned at frequency / =< / > /e, where < I > is the average current. If it possible 
to monitor these charge oscillations in time, one can actually realise an electron counter. 

The main problem for the realisation of such an electron counter is to develop a charge 
detector with a high charge sensitivity in combination with a large bandwidth. Recently a 
promising step forward has been made through the development of a radio frequency 
single-electron transistor (RF-SET). The working principle has already been described in 
bee. |L5]|. Possible error mechanisms in the electron counter based on RF-SET have yet 



to be analyzed. 



7 Conclusions 

We have reviewed the status of the understanding of single-electron devices in metrology. 
By pumping single electrons one can generate a current which is determined by the pump 
frequency and the fundamental electron charge. There are two possible routes for utiliza- 
tion of this idea. One route is to devise a high precision current standard based on the 
pumping principle. This would also allow for a comparison with a current measured on 
the basis of a Josephson voltage standard and a quantum Hall resistance standard. The 
direct comparison of the two currents would give a check of the internal consistency of the 
quantum standards. There are several technical difficulties in this, one of the main prob- 
lems being the small currents (presently a few picoamperes) generated by single-electron 
pumps. 

The second possible route is to realize a new capacitance standard, circumventing the 
cumbersome use of the calculable capacitor. The idea is to pump a known charge onto a 
capacitor and then to measure the resulting voltage. Note that this is a dc capacitance 
standard. In order to link its value to that of the calculable capacitor, which operates at 
1592 Hz, more knowledge is required about the frequency dependence of these standards. 

The state of the art 7-j unction single-electron pump has an intrinsic error per pumped 



electron of 1 part in 10^ [17], while generating a picoampere current. This approaches 
the level of metrological accuracy currently obtained for Josephson voltage standards and 
quantum Hall standards. However, the measured error rate is many orders of magnitude 
larger than that predicted by theory. 
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Although the quantum Hall effect is also not completely understood, it is an inherent 
quantum effect, whereas single-electron current generators are devices that have to be 
designed in order to suppress quantum mechanical effects (e.g. cotunneling) . Furthermore, 
the quantized Hall resistance has been shown to be independent of material parameters, 
while single-electron generators are very sensitive to the design layout and choice of device 
parameters. For these reasons a detailed understanding of the physical mechanisms of 
errors is essential for the use in metrology, and future research is needed in this direction. 
Therefore, at this stage, the single-electron pumps cannot be considered as being at a 
similar level as the Josephson voltage standard or the quantum Hall standard, neither 
experimentally nor theoretically. 

In this review, we have mainly concentrated on conventional tunnel junction SET de- 
vices such as pumps and turnstiles. Although the present accuracy is still not sufficient for 
all metrological applications, the ongoing technical developments justify the expectation 
of considerable improvements in device performance. For future designs several issues are 
important. For example to try waveforms which have been predicted to be much better 
than the triangular shape. Secondly, to minimize noise from charge traps in the substrate 
by alternative designs or materials. Also, the complexity of the design of turnstiles and 
pumps may be reduced by introducing high-ohmic microresistors in the bias leads thus 
suppressing cotunneling. In this way one could obtain an error per transferred electron 
equal to that of a 7-junction pump, but with less tunnel junctions. It makes the operation 
of the device easier and makes the device more robust. 

In addition, there are promising candidates for alternative single-electron current stan- 
dards. One example is the SAW device, another example is a single-electron counter based 
on a high-speed SET transistor. 

A large worldwide research is evolving in the field of manipulation of single electron 
charges. The efforts give enough confidence in a future realisation of a high precision SET 
current standard. 



8 Acknowledgments 



We are indebted to A. Zorin and D. Haviland for their input to Sees. 3.5 and 3.2 and to 



M. Keller for suggestions and careful reading of the manuscript. We thank J. Martinis, H. 
D. Jensen, and A. Hartland for discussions. This work was supported by the EU (SETamp 
project, contract number SMT4-CT96-2049). 



References 

[*] Present address: Niels Bohr Institute f. APG, 0rsted Laboratory, Universitetsparken 
5, DK-2100 Copenhagen 

[**] Permanent address: Nuclear Physics Institute, Moscow State University, Moscow 
119899 GSP, Russia. 

[1] H. R. Zeller and I. Giaver, Phys. Rev. 181, 789 (1969). 

[2] T. A. Fulton and G. J. Dolan, Phys. Rev. Lett. 59, 109 (1987). 

[3] D. V. Averin and K. K. Likharev, in Mesoscopic Phenomena in solids, edited by B. 
L.Al'tshuler, P. A. Lee, and R. A. Webb (Elsevier, Amsterdam, 1990). 



Towards single-electron metrology 



36 



[4] P. Delsing, K. K. Likharev, L. S. Kuzmin, and T. Claeson, Phys. Rev. Lett. 63, 1180 
(1989). 

[5] Y. V. Nazarov, Zh. Eksp. Teor. Fiz. 95, 975 (1989) [Sov. Phys. JETP 68, 561 (1989)]; 
Pis'ma Zh. Eksp. Teor. Fiz. 49, 1349 (1989) [JETP Lett. 49, 126 (1989)]. 

[6] M. H. Devoret et al, Phys. Rev. Lett. 64, 1824 (1990). 

[7] S. M. Girvin et al, Phys. Rev. Lett. 64, 3183 (1990). 

[8] K. Flensberg et al, Physica Scripta T42, 189 (1992). 

[9] D. B. Haviland, Y. Pashkin, and L. S. Kuzmin, Physica 203, 347 (1994). 

[10] L. S. Kuzmin et al, Phys. Rev. B 54, 10074 (1996). 



[11] W. Zheng et al, |cond-mat /9806044 . 



[12] P. Delsing, K. K. Likharev, L. S. Kuzmin, and T. Claeson, Phys. Rev. Lett. 63, 1861 
(1989). 

[13] L. J. Geerligs et al, Phys. Rev. Lett. 64, 2691 (1990). 
[14] H. Pothier et al, Europhys. Lett. 17, 249 (1992). 

[15] J. M. Martinis, M. Nahum, and H. D. Jensen, Phys. Rev. Lett. 72, 904 (1994). 

[16] M. W. Keller, J. M. Martinis, N. M. Zimmerman, and A. H. Steinbach, Appl. Phys. 
Lett. 69, 1804 (1996). 

[17] M. W. Keller, J. M. Martinis, A. H. Steinbach, and N. M. Zimmerman, IEEE Trans, 
lustrum. Meas. 46, 307 (1997). 

[18] M. W. Keller, J. M. Martinis, and R. L. Kautz, Phys. Rev. Lett. 80, 4530 (1998). 

[19] E. R. Wilhams, R. N. Ghosh, and J. M. Martinis, J. Res. Natl. Inst. Stand. Technol. 
97, 299 (1992). 

[20] K. K. Likharev and A. Zorin, J. Low. Temp. Phys 59, 347 (1985). 

[21] L. S. Kuzmin and D. B. Haviland, Phys. Rev. Lett. 67, 2860 (1991). 

[22] L. S. Kuzmin and D. B. Haviland, Physica Scripta T42, 171 (1992). 

[23] Y. Nazarov and A. A. Odintsov, in Single-electron tunneling and mesoscopic devices, 
edited by H. Koch and H. Liibbig (Springer- Verlag, Berlin, 1992), p. 91. 

[24] L. S. Kuzmin, Y. Pashkin, A. Zorin, and T. Claeson, Physica B 203, 376 (1994). 

[25] D. B. Haviland and P. Delsing, Phys. Rev. B 54, R6857 (1996). 

[26] J. M. Shilton et al, J. Phys.: Condens. Matter 8, L531 (1996). 

[27] V. I. Talyanskii et al, Phys. Rev. B 56, 15180 (1997). 

[28] N. M. Zimmerman, J. L. Cobb, and A. F. Clark, IEEE Trans. lustrum. Meas. 46, 
294 (1997). 



Towards single-electron metrology 



37 



[29] M. W. Keller, private communication. 

[30] B.N. Taylor and T.J. Witt, Metrologia 26, 47 (1989). 

[31] Single Charge Tunneling, edited by H. Grabert and M. Devoret, Nato ASI Series B 
Vol. 294 (Plenum Press, New York), 1992. 

[32] P. Lafarge et al, Z. Phys. B 85, 327 (1991). 

[33] T. A. Pulton, P. L. Gammel, and L. N. Dunkleberger, Phys. Rev. Lett. 67, 3148 
(1991). 

[34] A.F. Clark et al, Appl. Phys. Lett. 66, 2588 (1995). 
[35] M.J. Yoo et a/.. Science 276, 579 (1997). 

[36] A. N. Korotkov, M. R. Samuelsen and S. A. Vasenko, Journal of Applied Physics 76, 
3623 (1994). 

[37] C. Liu and Q. Niu, Physics Reports 286, 349 (1997). 
[38] M. Field et al, Phys. Rev. Lett 77, 350 (1996). 

[39] Y. Y. Wei, J. Weiss, K. von Klitzing, and K. Eberl, Appl. Phys. Lett. 71, 2514 (1997). 
[40] R. J. Schoelkopf et al. Science 280, 1238 (1998). 

[41] J. P. Pekola, K. P. Hirvi, J. P. Kauppinen, and M. A. Paalanen, Phys. Rev. Lett. 73, 
2903 (1994). 

[42] S. Farhangfar et al, J. Low Temp. Phys. 108, 191 (1997). 

[43] L. I. Glazman and K. A. Matveev, Pis'ma Zh. Eksp. Teor. Fiz. 48, 403 (1988) [JETP 
Lett. 48, 445 (1988).]. 

[44] D. V. Averin and A. N. Korotkov, zh. Eksp. Teor. Fiz. 97, 1661 (1990) [Sov. Phys. 
JETP 90, 937 (1990)]. 

[45] C. W. J. Beenakker, Phys. Rev. B 44, 1646 (1991). 

[46] V. F. Anderegg et al, Physica B 165-166, 61 (1990). 

[47] L. P. Kouwenhoven et al, Z. Phys. B - Condensed Matter 85, 367 (1991). 

[48] C. Urbina et al, IEEE Trans. Magn. 27, 2578 (1991). 

[49] S. Verbrugh, ph.D. thesis. Delft University of Technology. 

[50] H. D. Jensen and J. M. Martinis, Phys. Rev. B 46, 13407 (1992). 

[51] H. Pothier, Ph.D. thesis. University of Paris 6, 1991. 

[52] A. Iwasa, A. Fukushima, and A. A. Odintsov, Jpn. J. Appl. Phys. 34, 5871 (1995). 

[53] L. R. C. Fonseca, A. N. Korotkov, and K. K. Likharev, Appl. Phys. Lett. 69, 1858 
(1996). 



Towards single-electron metrology 



38 



[54] L. R. C. Fonseca, A. N. Korotkov, and K. K. Likharev, J. Appl. Phys. 79, 9155 
(1996). 

[55] D. V. Averin and A. A. Odintsov, Phys. Lett. A 140, 251 (1989). 

[56] D. V. Averin and Y. V. Nazarov, Phys. Rev. Lett. 65, 2446 (1990). 

[57] T. M. Eiles, G. Zimmerh, H. D. Jensen, and J. M. Martinis, Phys. Rev. Lett. 69, 148 
(1992). 

[58] C. Pasquier et al, Phys. Rev. Lett. 70, 69 (1993). 

[59] D. J. Paul, J. Cleaver, H. Ahmed, and T. WhaU, Phys. Rev. B 49, 16514 (1994). 
[60] H. Matsuoka and S. Kimura, Appl. Phys. Lett. 66, 613 (1995). 

[61] A. N. Korotkov, D. V. Averin, K. K. Likharev, and S. A. Vasenko, in Single- electron 
tunneling and mesoscopic devices, edited by H. Koch and H. Liibbig (Springer-Verlag, 
Berlin, 1992). 

[62] Y. V. Nazarov, J. Low Temp. Phys. 90, 77 (1993). 
[63] D. V. Averin, Physica B 194-196, 979 (1994). 
[64] K. Flensberg, Phys. Rev. B 55, 113118 (1997). 

[65] L. J. Geerligs and M. M. J. E. Mooij, Physica B 194-196, 1267 (1994). 

[66] J. Konig, H. Schoeller, and G. Schon, Phys. Rev. Lett. 78, 4482 (1997). 

[67] P. Joyez et al, Phys. Rev. Lett 79, 1349 (1997). 

[68] P. Lafarge and D. Esteve, Phys. Rev. B 49, 13309 (1993). 

[69] A.A. Odintsov, V. Bubanja, and G. Schon, Phys. Rev. B 46, 6875 (1992). 

[70] A. Zorin, private communication. 

[71] J. M. Martinis and M. Nahum, Phys. Rev. B 48, 18316 (1993). 

[72] J. M. Martinis, M. H. Devoret, and J. Clarke, Phys. Rev. B 35, 4682 (1987). 

[73] D. Vion et al., J. Appl. Phys. 77, 2519 (1995). 

[74] A. B. Zorin, Rev. Sci. Instrum. 66, 4296 (1995); D. C.Glatth et al., J. Appl. Phys. 
81, 7350 (1997). 

[75] V. Bouchiat, ph.D. thesis. University of Paris 6, 1997, Ch.5. 

[76] G. Zimmerh, T. M. Eiles, R. L. Kautz, and J. M. Martinis, Appl. Phys. Lett. 61, 237 
(1992). 

[77] A. B. Zorin et al, Phys. Rev. B 53, 13682 (1996). 
[78] E. H. Visscher et al, Appl. Phys. Lett. 66, 305 (1995). 



Towards single-electron metrology 



39 



[79] V. A. Krupenin, D. E. Presnov, M. N. Savvateev, H. Scherer, A. B. Zorin and 
J. Niemeyer, Conference on Precision Electromagnetic Measurements CPEM 1998, 
Conference Digist, ed. by T. L. Nelson, p. 140. 

[80] V. A. Krupenin, D. E. Presnov, M. N. Savvateev, H. Scherer, A. B. Zorin and 
J. Niemeyer, J. Appl. Phys, in press. 

[81] A. J. Manninen and J. P. Pekola, Czech. J. Phys. 46, 2293 (1996). 

[82] I. O. Kulik and R. I. Shehkter, Sov. Phys. JETP 41, 308 (1975). 

[83] For general information about the nanostructure simulators use http:/ /web- 



ext2.darpa.mil/ETO/ULTRA/Modeling.html . 



[84] R. H. Chen. MOSES: A general Monte Carlo simulator for single-electronic 
circuits. In Meeting Abstracts, volume 96-2, page 576. The Electrochem- 
ical Society, October 1996. For info about MOSES software package use 
ftp: I /hana.physics. sunysb.edu / pub / moses / . 

[85] L. R. C. Fonseca, A. N. Korotkov, K. K. Likharev, and A. A. Odintsov, J. Appl. 
Phys. 78, 3238 (1995). 

[86] C. Wasshuber, H. Kosina, Superlatt. Microstr., 21, No.l, 37 (1997); C. Wasshuber, 
About Single-Electron Devices and Circuits, Ph.D. Thesis, Institute for Microelec- 
tronics, University 



of Vienna, Austria, 1997, http://a63.iue.tuwien.ac.at/diss/wasshuber/diss/diss.html 



For info about SIMON software package use bttp: / /members. magnet. at / catsmeow/] . 

[87] N. S. Bakhvalov, G. S. Kazacha, K. K. Likharev, and S. I. Serdyukova. Sov. Phys. 
JETP, 68, 581-587 (1989); K. K. Likharev, N. S. Bakhvalov, G. S. Kazacha, and S. I. 
Serdyukova, IEEE Trans. Magn. 25, 1436 (1989). 

[88] M. Kirihara, N. Kuwamura, K. Taniguchi, and C. Hamaguchi, Extended Abstracts 
of the International Conference on Solid State Devices and Materials, pp. 328-330, 
Yokohama, Japan 1994. 

[89] S.A. Roy, Simulation Tools for the Analysis of Single Electronic Systems,PhD thesis. 
University of Glasgow, 1994. 

[90] D. V. Averin, A. A. Odintsov, and S. V. Vyshenskii, J. Appl. Phys. 73, 1297 (1993). 

[91] R. Bauerschmitt and Y. V. Nazarov, Phys. Rev. B 47, 9997 (1993). 

[92] J. M. Shilton et al, J. Phys.: Condens. Matter 8, L337 (1996). 

[93] L. J. Geerligs, et al, Z. Phys. B 85, 349 (1991). 

[94] A. A. Odintsov, Sov. Phys. JETP Lett. 60, 738 (1994). 



